SOME CONSIDERATIONS OF MATRIX EQUATIONS 
USING THE CONCEPT OF REPRODUCTIVITY 

BRANKO MALESEVICi AND BILJANA RADICIC^ 

Abstract. In this paper we analyse Cline's matrix equation, generalized Pen- 
rose's matrix system and a matrix system for /c-commutative {l}-inverses. We 
determine reproductive and non-reproductive general solutions of analysed matrix 
equation and analysed matrix systems. 



1. Introduction 

In this paper we determine general and reproductive general solutions of analysed 
matrix equation and analysed matrix systems. We are going to use the concept of re- 
productivity in order to prove that certain formulas represent the general solutions of 
analysed matrix equation and analysed matrix systems. The concept of reproductive 
equations was introduced by S.B. Presic [2] in 1968. 

Let iS be a given non-empty set and J be a given unary relation of S. Then an equa- 
tion J{x) is consistent if there is at least one element xq G S, so-called the solution, 
such that J{xq) is true. A formula x = 4>{t), where : 5 — )■ is a given function, 
represents the general solution [20J of the equation J{x) if and only if 

(Vt)J((/)(t)) A (Vx)(J(x) =^ (3t)x = 4){t)). 

Let us cite the definition of reproductive equations according by S.B. Presic [2]. 

Definition 1.1. The reproductive equations are the equations of the following form: 

X = ^p{x), 

where x is a unknown, S* is a given set and if : S — > is a given function which 
satisfies the following condition: 

(1.1) (po(f = if. 
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The condition (11. ip is called the condition of reproductivity [2J. The fundamental 
properties of reproductive equations are given by the following two statements (S.B. 
Presic [2]) (see also [5], [6] and [H]). 

Theorem 1.1. For any consistent equation J{x) there is an equation of the form 
X = if{x), which is equivalent to J{x) being in the same time reproductive as well. 

Theorem 1.2. If a certain equation J{x) is equivalent to the reproductive one x = f{x), 
the general solution is given by the formula x = if{y), for any value y & S. 

Let us remark that a formula x = 4>{t), where (p : S ^ S is a given function, represents 
the reproductive general solution [20] of the equation J{x) if and only if 

(Vt)J(0(t))A(Vt)(J(t)^t = 0(t)). 

Reproductivity of some equations of mathematical analysis was studied by J.D. 
Keckic in [9], [To]. In [15] J.D. Keckic and S.B. Presic considered the gen- 
eral applications of the concept of reproductivity. The general applications of the 
concept of reproductivity in various mathematical structures can also be found in 

[i, m and m- 

2. Main results 

Let m,n E N and C is the field of complex numbers. The set of all m x matrices 
over C is denoted by C*"^". By C^^" we denote the set of all matrices from C™^"" 
with a rank a. For AgC"^" the rank of A is denoted by rank{A). The unit matrix of 
order m is denoted by (if the dimension of unit matrix is known from the context, 
we omit the index which indicates the dimension and we use designation /). Let 
A G C™^", then a solution of the matrix equation 

AX A = A 

is called {l}-inverse of A and it is denoted by A^^K In the general case {l}-inverse 
of A is not uniquely determined. The set of all {l}-inverses of A is denoted by A{1}. 
It can be shown that ^4(1} is not empty. {l}-inverse of A is uniquely determined if 
A is regular. In that case {l}-inverse A^^^ corresponds to A~^ i.e. ^4(1} = 
There are also other types of inverses. More informations about {l}-inverse and 
other types of inverses can be found in [18] and [T9j. For A G C™^"^ the smallest 
non-negative integer k such that rank{A^) = rank{A^^^) is called the index of A and 
it is denoted by Ind{A). 
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This section of paper is divided into three parts. The first part is devoted to the 
matrix equation 

(2.1) A^XB'' = C, 

where A e C^^p, B e C G C^^", m > k = ind{A) and n > / = ind{B). In the 

second part we consider the matrix system 

(2.2) (2.2.1.) A'^X = B A (2.2.2.) XD"" = E, 

where A e C^^^, B e C^^^, D e C'^'"', E e Cp^«, m > Ind{A) and n > Ind{D). A 
solution of the matrix system 

(2.3) AXA = A A A''X = XA'', 

where A G C^^p is a singular matrix and A; G N, is analysed in the third part of this 
section. 

2.1. In this part we analyse the matrix equation (12.11) . In the paper [3] R.E. Cline 
was the first one who considered the matrix equation (12. ip . Using Penrose's con- 
dition for the consistence of the matrix equation AXB = C, R.E. Cline concluded 
that the matrix equation (12. ip is consistent if and only if 

(2.4) A™(A'")WC(S")WS" = C. 

In the paper [3] it was shown that the matrix equation (12.11) is consistent for any 
m > k and any n > / if and only if the matrix equation A^XB^ = C is consistent. 
Based on the results in the paper [22] the condition of consistence (12. 4p for the matrix 
equation (12. ip can be also considered in a new form (see Theorem 2.1). 

Lemma 2.1. If the matrix equation (12. ip is consistent, the equivalence 

^m^Qu = c ^ X = f{X) = X - (A'")(i)(A'"X5" - C){B''Y^^ 

is true. 

Proof. : Suppose that A^XB^ = C. Then, the equality 

(A")(i)A"XS"(S")(^) = (A")(i)C(5")(^) 

is also true and 
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= X- (A'")(i)(A™X5" - C)(5")(i) 

= /W 

=) : Suppose that X = f{X) = X - {A'^Y'^\A^ X B"' - C){B"Y^\ Then, 
= A"'f{X)B'' 

^ V ' V ' ^ V ' 

{=A-) (=S") (^C) 

= - ^'"XE" + C 

= C. 



□ 



Remark 2.1. It is easy to show that piY) = f{Y) i.e. the function / satisfies the 
condition of reproductivity. Therefore, if the matrix equation (12.11) is consistent, it is 
equivalent to the reproductive matrix equation X = f{X). 

Based on the previous remark and Theorem 11.21 we conclude that the following theo- 
rem is true. 

Theorem 2.1. If the matrix equation (12. ip is consistent, the general solution of the 
matrix equation (12. ip is given by the formula 

where Y is an arbitrary matrix corresponding dimensions. 
The following theorem is an extension of the previous theorem. 

Theorem 2.2. If Xq is a particular solution of the matrix equation (12.11) . the general 
solution of the matrix equation (12.11) is given by the formula 



(2.5) X = g{Y) = Xo + Y- (A™)(i)A™yS"(S")(i), 

where Y is an arbitrary matrix corresponding dimensions. 
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Proof. It is easy to see that the solution of the matrix equation (12.11) is given by (12. 5p . 
On the contrary, let X be any solution of the matrix equation (12. ip . then 

= X- + (yl™)(^)A™Xo5"(S")(^) 

= X- (A'")(iM™(X - Xo)B"(S")(^) 

= Xo + {X- Xq) - {A'''Y^'>A"'{X - Xo)B''{B'''f^ 

= Xo + Y- 

= 9iY), 

where Y = X — Xq. From this we see that every solution X of the matrix equation 
( 12. ip can be represented in the form (12. Sp . □ 

Remark 2.2. From g'^{Y) = g{Y) + [Xq - (A"")(i)C(S")(^)) we conclude that the 
function g is reproductive if and only if Xq = (A"^)^^^C{B"')^^K 

Remark 2.3. Theorem 12.21 is an extension, as we mentioned, of Theorem 12.11 be- 
cause there is a matrix equation (12. ip and a particular solution Xq such that Xq ^ 
(^m)(i)(7(^")(i) for any choice of {l}-inverses (A™-)*^^) and (i?")^-'^) similar to the cor- 
responding example from [21] . 



2.2. In this part we analyse the matrix system (12. 2p . as a special extension of 
Penrose's matrix system [1]: 

AX = B A XD = E, 

using the concept of reproductivity. 

Based on the result from we conclude that one common solution of the matrix 
system (12.20 is given by 

= + ^(Z}")(i) - (A™)(i)(A")A^(D")(i). 

The results which follow are extensions of the results from [18] (pp. 54-55) and [22] . 

Lemma 2.2. The matrix equations (2.2.1.) and (2.2.2.) have a common solution if 
and only if each equation separately has a solution and 

A'^E = BD"". 
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Proof. The proof is similar to the proof in |18] • 



□ 



Lemma 2.3. If the matrix system (12. 2p is consistent, the equivalence 
(A^X = B A XD" = E) 

X = f{X) = Xi + (/ - - 



is true. 



Proof. The proof is similar to the proof in [22j . 



□ 



Remark 2.4. It is easy to show that f^{Y) = f{Y) i.e. the function / satisfies the 
condition of reproductivity. Therefore, if the matrix system (12. 2p is consistent, it is 
equivalent to the reproductive matrix equation X = f{X). 

Based on the previous remark and Theorem 11.21 we conclude that the following theo- 
rem is true. 

Theorem 2.3. // the matrix system (12.21) is consistent, the general solution of the 
matrix system (12.21) is given by the formula 



where Y is an arbitrary matrix corresponding dimensions. 

The following theorem is an extension of the previous theorem. 

Theorem 2.4. If Xo is a particular solution of the matrix system (12.21) . the general 
solution of the matrix system (12. 2p ) is given by the formula 



X = f{Y) = Xi + (/ - - 



X = g{Y) =Xo + {I- (A'")(i)A'")y(/ - 



where Y is an arbitrary matrix corresponding dimensions. 



Proof. The proof is similar to the proof in P2]- 



□ 



Remark 2.5. From g'^{Y) = g{Y) + {Xq — Xi) we conclude that the function g is 
reproductive if and only if Xq = Xi. 
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2.3. In this part we analyse the matrix system (12. 3p . The second equation of the 
matrix system fl2.3p determines {5^}-inverse of A. A solution of the matrix system 
(12. 3 p is {1, 5'^}-inverse which is called k- commutative {l}-inverse and is denoted by A. 
fc-commutative {l}-inverses were considered in [11], [12] and [16]. It is easy to check 
that one solution of the matrix system (12. 3p is given hy X = AAA. In [12J J.D. 
Keckic gave the condition for the consistency of the matrix system (12. 3p . We are 
going to represent the formula of the general reproductive solution for the consis- 
tent matrix system (12. 3 p using the concept of reproductive equations. We need the 
following four lemmas. 

Lemma 2.4. A^A^ = A^A^. 

Proof. A^A^= A^AA^-^=AA''A^-^ = AA^A''-'^= A^A^A^-^ = ... =A''A''. □ 

(=AA*=) (=AA'') 

Lemma 2.5. A'^A'^A'' = A'' . 

Proof. Let us note A'' = A'=~M^= A''-^AAA = A^AA = AA''+\ Therefore, 

{=AAA) (=1a'=) 

A^ = A{A'') A = A{AA''+^) a = A^A^A^ = ... = A^A^A^ ^ ^ A'' A'' A''. □ 

Lemma 2.6. For any particular solution Xq of the matrix system (12. 3 p the equalities 
XqA^A^ = A^A^+^ and A^A^Xq = A^A^^^ 

are true. 

Proof. We are going to prove the first equality. 

Xq ^ A^ = XqA'' A^ A'' = A^XqA^A"^^ 

= A''-^ AXoAA^-^A^^ = A^-^AA^-^A^^ = A^^-^A^^, 

{=A) 

j^k — j\k A^ A^ A^^'^ = y42A;^2fc+l _ j\k AA^^ = A'^ AA^ A^'^ 

( = A'^A'^A'') ( = {=AA'') 

= A^'^ AAAA^~^A^^ = A^~^ AA^~'^ A^^ = 

i=A) 

Hence, XqA^A'^ = A'^A''^^. The second equality is proved similarly. □ 
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Lemma 2.7. Let X = AAA. If the matrix system ( \2.3h is consistent, the equivalence 
AX A = A A A^X = XA'' 

X = f{X)=X+X-{I-AA)XA''A''-A''A''X{I-AA)- AaxaA 

is true. 

Proof. ^) : Suppose that AX A = A A A'^X = XA\ Then, 

AAXA'^A'' = AAXAA^-^A^ = A^A^ = A'^AA^ = A''A''+\ 

(=A) {=yl*A) 

Bearing in mind that XA^A^ = A'^A'^^^ (Lemma I2.6p . we conclude that 

XA'^A'' = AAXA^A^. 

In a similar way we get that A^A^X = A^A^XAA. 
Therefore, 

X = AAA + X - XA'^A'' + AAXA'^A^ - A'^A'^X + A'^A'^XAA - A^A 

(=AXA) 

= AAA + X - XA^A^ + AAXA^A^ - A^A^X + A^A^XAA - AAXAA 
= AAA + X -{I- AA)XA^A^ - A^A^X{I - AA) - AAXAA 
= fiX). 

^) : Suppose that X = f{X) = AAA + X - {I - AA)XA^A^ - A''A^X{I - AA) 
-AAXAA. Then, 

AX A = Af{X)A 

= A{AAA + X -{I- AA)XA''A'' - A^A''X{I - AA) - AAXAA) A 

= AAAAA+AXA-A{I-AA) XA^A^A-AA^A^X {I-AA)A-AAAXAAA 

i=A) ' J^) ' ' J^) ' {=A) (=A) 

= ^AA^+AXA - AX A 

(=A) 
= A, 
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A'^X = Af^fiX) 

= A'' {AAA + X - (/ - AA)XA''A'' - A^A^X{I - AA) - AAXAA) 

= A'' AAA + A''X-A\I~ AA)XA^A'' -^A^J^X {I-AA)-A''A AXA A 

{ = Ak) {=A) 

YES] 

= A^AAA + A^X-A'^XA^A^ + A^ AAX A'' - A'' X + A'' X AA- A'' AaA 

= -A'' X/^A^ +a^Aaxaa^-^A^ + a^-^axaA 

( = AfeAfe+i) {=A) {=A) 

LEU 

= -A^A^A^+^ + A'^AA^A^ + A^A 

{=A*'A) 

= -A2fc^fc+i _^ A^Af'AA^ + A^A 

= A'' A, 
XA^ = f{X)A^ 

= {AAA + X-{I- AA)XA^A^ - A^A''X{I - AA) - AAXAA)A^ 

= AAAA'' +XA^-{I-AA)X A^A^-A'' A''X{I~AA)A''-A AXA AA^ 

( = A*) (=A) 
^L[23] ' 

= AaAA^ +XA^-XA^+ a ax a'' - A'^A^XA^ + A^A^XAAA^ - AaAA^ 
= AAXAA'''^ - A^A^ A'' + A^A^XAAA^ 

{=A) ( = A'=A'=+i) 

= AA'' - A'^A'^+^A'' + A^A^XAAA^ 
= AA^ - A^A^+^A^ + A^A^-^ AXAAA^ 

i=A) 

= AA'' - A'^A'^+^A'' + ^A^ AA'' 

( = AfcAfe) 

= AA'' — A'^A'^^^A'' + A'^A'^^'^A'^ 
= AA''. 

From A'' A = AA'' we see that A'^X = XA''. □ 



10 



BRANKO MALESEVIC AND BILJANA RADICIC 



Remark 2.6. It is easy to show that f^(Y) = f{Y) i.e. the function / satisfies the 
condition of reproductivity. Therefore, if the matrix system (12. 3p is consistent, it is 
equivalent to the reproductive matrix equation X = f{X). 



Based on the previous remark and Theorem 11.21 we conclude that the following theo- 
rem is true. 

Theorem 2.5. // the matrix system (12.30 is consistent, the general solution of the 
matrix system (12. Sp is given by the formula 

X = f{Y) =AAA + Y - {I-AA)YA^A^ - A'^A^Yi.I-AA) - AAYAA, 

where Y is an arbitrary matrix corresponding dimensions. 

In [12] J.D. Keckic also proved this theorem, but his proof is different from the 
previously exposed proof. 

The following theorem is an extension of the previous theorem. 

Theorem 2.6. If Xq is a particular solution of the matrix equation (12.31) . the general 
solution of the matrix equation (12.31) is given by the formula 



(2.6) X = g{Y) =Xo + Y - {I - AA)YA''A'' - A'^A'^Y {I - AA) - AAYAA. 
where Y is an arbitrary matrix corresponding dimensions. 

Proof. It is easy to see that the solution of the matrix system (12. 3p is given by (12. 6p . 
On the contrary, let X be any solution of the matrix equation (12. 3p . then 



X = X- A^A^ + A^A^ +^^'=+^'=+1 - AA''+^A''+^ - A^A^ 

( = XAkAk) ( = XqA'^A'') ( = ^feyl'=J 

L iZGl LiMl LED 

+ AA + A'^A''^^ - A''A''+^ AA-AAA + AAA 



( = AkA^^X) ( = A'^A'^Xq) ( = Afe^^Xo) {=AXA) {=AXqA) 



■ X - XA'^A^ + XqA^A^ + AA 4^4^ - AA A^4^ - A^A^X 

( = XAi'Ai') ( = XoA^yifc) 

+ A^A^XAA + A^A^Xq - A^A^XqAA - AAXAA + AAXqAA 
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= Xo+iX-Xo)-{X-Xo)A''A^+AAXA''A''-AAXoA^A''-A^A\X-Xo) 

+ A^A'^i^X - Xo)AA - AA{X - Xq)AA 
= Xq + {X- Xo) -{X- Xq)A''A'' + AA{X - Xq)A^A^ - A'^A'^i^X - Xo) 

+ A^'A'^i^X - Xq)AA - AA{X - Xq)AA 
= Xo + (X - Xo) -{I- AA){X - Xq)A^A^ - A^A^{X - Xo)(/ - AA) 

- AA{X - Xa)AA 
= Xq + Y - {I-AA)YA^A^ - A^A^Y{I-AA) - AAYAA 
= 9iY), 

where Y = X — Xq. From this we see that every solution X of the matrix system 
( 12 ■2p can be represented in the form (I2.6p . □ 

Remark 2.7. From g'^{Y) = g{Y) + {Xq — X) we conclude that the function g is 
reproductive if and only if Xq = X. 

Remark 2.8. The preceding result is an extension of the consideration which is given 
in [22] (see Application 2.2) 

3. Conclusion 

We want to emphasize that there are also other matrix equations and matrix systems 
whose solutions can be analysed in the same way as we done with fl2.ip . fl2.2p and 

(I23D. 

Acknowledgment. The research is partially supported by the Ministry of Science 
and Education, Serbia, Grant No. 174032. 

References 

[1] R. Penrose, A generalized inverses for matrices, Math. Proc. Cambridge Philos. Soc. 51 (1955), 
406-413. 

[2] S.B. Presic, Une classe d'equations matricielles et I'equation fonctionnelle — /, Pub- 
lications de I'institut mathematique, Nouvelle serie, tome 8 22, Beograd 1968, 143-148. 
(http://publications.mi.sanu.ac.rs/) 

[3] R.E. Cline, Inverses of rank invariant powers of a matrix, SIAM J. Numer. Anal., 5 (1968), 
182-197. 

[4] S.B. Presic, Une methode de resolution des eequations dont toutes les solutions appartiennent 
d un ensemble fini donne, C. R. Acad. Sci. Paris Ser. A 272 (1971) 654-657. 



12 



BRANKO MALESEVIC AND BILJANA RADICIC 



[5] S.B. Presic, Ein Satz Uber Reproduktive Losungen, Publications de I'institut mathematique, 
Nouvelle serie, tome 14 28, Beograd 1972, 133-136. 

[6] M. Bozic, A Note On Reproductive Solutions, Publications de I'institut mathematique, Nouvelle 
serie, tome 19 33, Beograd 1975, 33-35. 

[7] S. Rudeanu, On reproductive solutions of arbitrary equations, Publications de I'institut mathe- 
matique, Nouvelle scric, tome 24 38, Beograd 1978, 143-145. 

[8] D. Bankovic, On general and reproductive solutions of arbitrary equations, Publications de 
I'institut mathematique, Nouvelle serie, tome 26 40, Beograd 1979, 31-33. 

[9] J.D. Kcckic, Reproductivity of som,e equations of analysis I, Publications de I'institut mathe- 
matique, Nouvelle serie, tome 31 45, Beograd 1982, 73-81. 
[10] J.D. Keckic, Reproductivity of some equations of analysis II, Publications de I'institut mathe- 
matique, Nouvelle serie, tome 33 47, Beograd 1983, 109-118. 
[11] L.D. Dobryakov, Commuting generalized inverse matrices. Mathematical Notes, Volume 36, 
Number 1, 500-504, 1985. (Translated from Matematicheskie Zametki, Vol. 36, No. 52, 17-23, 
1984.) 

[12] J.D. Keckic, Commutative weak generalized inverses of a square matrix and some related matrix 
equations. Publications de I'institut mathematique, Nouvelle serie, tome 38 52, Beograd 1985, 

39-44. 

[13] D. Bankovic, All solutions of finite equations. Discrete Mathematics, Volume 137, 1995, 1-6. 
[14] D. Bankovic, Formulas of general solutions of Boolean equations. Discrete Mathematics, Volume 
152, 1996, 25-32. 

[15] J.D. Keckic and S.B. Presic, Reproductivity - A general approach to equations. Facta Universi- 
tatis (Ms), Ser. Math. Inform. 12 (1997), 157-184. 

[16] J.D. Keckic, Some remarks on possible generalized, inverses in semigroups. Publications de 
I'institut mathematique, Nouvelle serie, tome 61 75, Beograd 1997, 33-40. 

[17] S.B. Presic, A generalization of the notion of reproductivity. Publications de I'institut mathe- 
matique, Nouvelle serie, tome 67 81, Beograd 2000, 76-84. 

[18] A. Ben-Israel and T.N.E. Greville, Generalized Inverses: Theory and Applications, Springer, 
2003. 

[19] S.L. Campbell and CD. Meyer, Generalized Inverses of Linear Transformations, Society for 

Industrial and Applied Mathematics, 2009. 
[20] D. Bankovic, General Solutions of System of Finite Equations, Scientific Publications of the 

State University of Novi Pazar Ser. A: Appl. Math. Inform, and Mech. vol. 3, 2 (2011), 117- 

121. 

[21] B. Malesevic and B. Radicic, Non-reproductive and reproductive solutions of some matrix equa- 
tions. Proceedings of the Conference Mathematical and Informational Technologies, MIT-2011, 
V. Banja, Serbia, 2011, 246-251. (http://mit.rs/, http://conf.nsc.ru/MIT-2011/) 

[22] B. Radicic and B. Malesevic, Some considerations in relation to the matrix equation AXB = C 
(http://arxiv.org/abs/1108.2485) 



^Department of Applied Mathematics 

Faculty of Electrical Engineering, University of Belgrade, 

Serbia 

E-mail address: malesevic@etf.rs 
^ part-time job 

Department of Mathematics, Physics and Descriptive Geometry 
Faculty of Civil Engineering, University of Belgrade, 
Serbia 

E-mail address: radicic. biljanaQyalioo. com 



